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^ ! Abstract 

W)' The heats of formation of HCIO4 and CI2O7 have been determined to chemical accuracy for 

' the first time by means of Wl and W2 theory. These molecules exhibit particularly severe 

degrees of inner polarization, and as such obtaining a basis-set limit SCF component to the 
total atomization energy becomes a challenge. (Adding high-exponent d functions to a standard 
spd basis set has an effect on the order of 100 kcal/mol for CI2O7.) Wilson's aug-cc-pV(n+d)Z 



O 



basis sets represent a dramatic improvement over the standard aug-cc-pVnZ basis sets, while 
the aug-cc-pVnZ+2dlf sequence converges still more rapidly. Jensen's polarization consistent 
basis sets still require additional high-exponent d functions: for smooth convergence we suggest 
the {aug-pcl+3d,aug-pc2+2d,aug-pc3+d,aug-pc4} sequence. The role of the tight d functions 
^ • is shown to be an improved description of the CI (3d) Rydberg orbital, enhancing its ability 

to receive back-bonding from the oxygen lone pairs. In problematic cases like this (or indeed in 
' general), a single SCF/aug-cc-pV6Z+2dlf calculation may be preferable over empirically motivated 

■ extrapolations. Our best estimate heats of formation are /S.H°^ 2m\^^^^i{^] — — O.Gitl kcal/mol 

f-*) , and /S.H°j 2Q^[C\20-j{g)\ = 65.9±2 kcal/mol, the largest source of uncertainty being our inability to 

\ account for post-CCSD(T) correlation effects. While G2 and G3 theory have fairly large errors, 

Q ' G3X theory reproduces both values to within 2 kcal/mol. 
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I. INTRODUCTION 



The author's Dirac medal lecture at WATOC'05 focused on some recent advances in t 



area o 

n 



computational thermochemistry in general, and on the Wn (Weizmann-n) theory 



ne 



, y developed in our own laboratory in particular. The subject has been reviewed very 
and therefore will not be re-reviewed here. Rather, we will present an example 
practical application. 

Perchloric acid, HCIO4, and its anhydride, CI2O7, have been known for nearly two 



centuries p]. Industrial uses of perchloric acid are manifold, and higher ch 



orine oxides 
(See also 



01 



have recently been inmlicated in theories of stratospheric destruction of ozone 
two recent papers 0, on the Cl20„ and Br20„ series (n=l-4) and their anions.) Yet no 
reliable thermochemical data are available: the reason for this is probably best illustrated 



by the following quote from a spectroscopic study 10| : 



Caution! Several explosions occurred during the course of this work. It was 
necessary to perform all experiments wearing heavy gloves and a face shield with 
the sample properly shielded. 



Cioslowski et al. compiled | 111] an extensive and chemically diverse set of thermochemical 
data for the purpose of benchmarking, calibration, and parametrization of computational 
thermochemistry methods. Their dataset includes CI2O7, for which an unsourced value 
of AiJj298~65.0 kcal/mol was taken from the Wagman et al. compilation 1^ . An atom- 
equivalent scheme based on B3LYP/6-311++G** calculations proposed in Ref. yielded 
its single largest error (91 kcal/mol) for this molecule. While it is to be expected that this 



it must also be said that 



molecule will exhibit very strong inner polarization effects 
the 65.0 kcal/mol number is a crude estimate, and that any comparison with it will be 
semiquantitative at best. 

For Ha04, Francisco Q ( motes Aifj298~4±4 kcal/mol from the NIST database^^, and 
a group additivity estimate|l7l| of -1.5 kcal/mol. Francisco himself calculated AifjQ=-|-10.8 
kcal/mol using G2 theory. Similar remarks apply as for CI2O7. 

In the present work, we will apply Wl and W2 theor y to these systems. We will also show 
that they are extreme cases of inner polarization effectsjisl. Q] and that any calculation that 
does not adequately take this into account is doomed to failure. 



2 



II. COMPUTATIONAL DETAILS 



All calculations were run on a 4-CPU AMD Opteron 846 server with 8 GB of RAM, 
two SATA system disks, and eight 72 GB UltraSCSI320 scratch disks, running SuSE Linux 
Enterprise Server 9 and custom-built for our group by Access Technologies of Rehovot, Israel. 
The scratch disks are operated as two hardware RAID-0 arrays of four disks aggregated in 
software using the Linux "md" facility. For large soft-RAID "chunk sizes" of 512 KB and 
up, we are able to obtain sustained streaming read and write bandwidth in excess of 300 
MB/s, as measured by IOzone |l8| for a 16 GB file. 

The CCSdQ and CCSD(T)0,|2l| coupled cluster calculations involved in Wl and W2 
theory were carried out using MOLPRO 2002. 6j2^ (always using conventional rather than 
integral-direct algorithms), while all remaining calculations were carried out using a locally 
modified version of Gaussian 03 2^. NBO (natural bond orbital) 24 1 and AIM (atoms- in- 
molecules) j^l analyses were carried out using the relevant modules of Gaussian 03. 

Three families of basis sets were used. The first are the original aug-cc-pVnZ (augmented 
correlation consistent polarized n-tuple zeta) basis sets of Dunning and coworkers[3]: 
besides the unmodified aug-cc-pVnZ sequence, we considered the aug-cc-pVnZ+2dlf 
sequence proposed in the original W1/W2 paper [l|, where the"+2dlf" notation represents 
the addition of two additional d and one additional / function, with exponents obtained by 
successively multiplying the highest exponent of that angular momentum already present 
by a factor of 2.5. 

The second family are the aug-cc-pV(n+d)Z basis sets of Wilson, Peterson, and 



Dunning 



27l |. which were specifically developed to cope with inner polarization effects and 



contain an additional d function beyond the original aug-cc-pVnZ basis sets. 



The third fami 
Jensen 



29 
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y are the aug-pcn (augmented polarization consistent) basis sets of 
. Unlike the correlation consistent basis sets which were developed 



for optimum recovery of the valence correlation energy, the polarization consistent basis 
sets were developed for optimum recovery of Hartree-Fock and DFT energies, for which 
requirements (in terms of optimum exponents and choice of polarization functions) are 
considerably different. We considered the aug-pcl, aug-pc2, aug-pc3, and aug-pc4 basis sets, 
which are of 3s2p, 4s3p2(i, 6s5p3d2f, and 8s7p4:d3f2g quality, respectively, on hydrogen, of 
4s3p2d, 5s4p3d2f, 7s6p5d3f2g, and 9sSp7d4:f3g2h quality, respectively, on oxygen, and of 
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5s4:p2d, 6s5p3d2f, 7s6p5d3f2g, and 8s7p7d4:f3g2h quality, respectively, on chlorine. Here 
too, we considered addition of one or more high-exponent d and / functions (denoted by 
"-|-d", "+2d" suffixes), with exponents again in even-tempered sequences with a stride 
factor of 2.5. 

Detailed descriptions (and rationales) for the different steps in Wl and W2 theory can be 
found in the original references 0,0,3,0]. In the interest of making the paper self-contained, 
we will briefly summarize them: 



a B3LYP/cc-pV(T+d)Zj22] reference geometry is used for Wl theory. Normally 
a CCSD(T)/cc-pV(Q+d)Z reference geometry is called for in W2 theory: as full 
CCSD(T)/cc-pV(Q+d)Z optimizations proved computationally intractable, we have 



instead used the best available DFT geometries 



3 (at the B97-l/aug-pc3+dlevei0)- 



The actual geometries used can be found in Figure 1; 

the SCF component is extrapolated from aug-cc-pVTZ-|-2dlf and aug-cc-pVQZ+2dlf 
basis sets in Wl theory, and from aug-cc-pVQZ+2dlf and aug-cc-pV5Z+2dlf basis 
sets in W2 theory, using the formulaj^ E{L) = E^o + A/L^; The acronyms Wlw and 
W2w (rather than Wl and W2) indicate that aug-cc-pV(n+d)Z rather than aug-cc- 
pVnZ+2dlf basis sets were used0|; 

the CCSD valence correlation component is extrapolated from the same pairs of basis 
sets using E(L) = E^r, + A./L°', where a = 3.22 for Wl and Wlw theory and a = 3 



for W2 theory [1 



the (T) valence correlation component is extrapolated using the same expression, but 
from aug-cc-pVDZ+2d and aug-cc-pVTZ+2dlf basis sets in the case of Wl theory, 
and aug-cc-pVTZ+2dlf and aug-cc-pVQZ-|-2dlf basis sets for W2 theory (with the 
appropriate substitutions in the case of Wlw and W2w). Note that the expensive (T) 
step need not be carried out in the largest basis set used overall. Note also that ROHF 



references are used for the open-shell atoms, and that the definition of Ref. 2l| for the 
ROHF-CCSD(T) energy is used; 



the differential contribution of inner-shell correlation is obtainec 
with a specially developed MTsmall core correlation basis set 



at the CCSD(T) level 

0Q; 
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the scalar relativistic contribution can be obtained either as first-order mass-velocity 



and Darwin corrections 



37l | at the ACPF (averaged coupled pair functional |38|) level 
with the MTsmall basis set, or (in the Gaussian 03 23] implementation of Wl theory 
and in the present work) as the difference between nonrelativistic and second-order 
Douglas-Kroll-HessQ CCSD(T) /MTsmall energies; 

These being closed-shell molecules, first-order spin-orbit corrections only affect the 
separated atoms and are obtained from the well-established atomic fine structures; 

The zero-point vibrational energy was obtained from scaled B3LYP/cc-pV(T+d)Z 
harmonic frequencies for Wlw theory, and from a recent DFT anharmonic force field 
study for W2w theory; 

Thermal corrections were obtained from the same data, although corrections are 
required for the hindered rotations in both molecules (see below). 



III. RESULTS AND DISCUSSION 



The reference geometries used are given in Figure 1. Most salient results can be found in 
Table 1: basis set convergence for the SCF component is presented in detail in Table 2. We 
shall discuss these results component by component. 



A. Hartree-Fock component 

Normally, the Hartree-Fock component is the easiest one to get right in an ab initio 
thermochemistry calculation. Yet in this case, very strong inner polarization effects cause a 
sharp dependence of the Hartree-Fock component of the binding energy on the presence of 
high-exponent d functions in the basis set. This can be seen in more detail in Table 2. 

For HCIO4 with the aug-cc-pVDZ or aug-pcl basis sets, addition of four high-exponent d 
functions (with exponents in an even-tempered series with stride 2.5) increases the computed 
binding energy by over 50 kcal/mol: for CI2O7, the magnitude of the effect reaches a 
whopping 100 kcal/mol! 

With aug-cc-pVTZ (or aug-pc2) and spdfg basis sets, the contributions for CI2O7 go down 
to 45 and 23 kcal/mol, respectively — the aug-pc3 basis set already contains one additional 
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high-exponent d function to begin with, and thus requires less additional d functions to 
reach saturation. 

An NBO analysis of the wave function reveals that, while the chlorine d orbitals have 
natural populations in the 0.35 range, these orbitals have Rydberg character and do not 
participate in any natural bond orbitals: as shown repeatedly previously 0, |^, there is 
no basis for describing HCIOa, nor by extension CI2O7, as hypervalent molecules. NPA, 

rr 

APT (Atomic Polar Tensor [42]), and AIM all yield partial charges consistent with a general 
H0(Cl+^)(-0^)3 bonding picture (see below). 

What happens to the wavefunction when extra d functions are added? Comparing NBO 
analyses for HCIO4 with aug-pcl and aug-pcl+Sd basis sets, we see that the energy of the 
3d Rydberg orbital on CI is lowered from 0.985 to 0.839 a.u. for dxy, 0.974 to 0.830 a.u. 
for d^z, 1-200 to 1.016 a.u. for dy^, 1.134 to 0.966 a.u. for dx2_y2, and 1.173 to 0.996 a.u. 
for (i322_^2. Thus the orbitals become more accessible: overall 3d population increases from 
0.28 to 0.37 electron equivalents, with the corresponding NBO occupations clustering in a 
"quasi-tas" triad of 0.117, 0.103, and 0.100 and a "quasi-Cg" dyad of 0.0594 and 0.0588. 
(These numbers add up to more than 0.37 as there are small non-d components in these 
orbitals.) Interestingly, the Wiberg bond indices for the three short Cl-0 bonds go up by 
about 0.05, and of the long bond by 0.02. The increase in the AIM (Bader) covalent bond 
orders is more pronounced: from 1.27 to 1.42-1.43 for the short bonds, and from 1.00 to 1.05 
for the long bond. Second-order perturbation theory analysis of the NBOs reveals, besides 
the expected interactions between oxygen lone pairs and antibonding CIO orbitals, strong 
interactions between oxygen lone pairs and the 3d Rydberg orbitals on chlorine. 

Consistent with all of this, the AIM charge distribution is considerably less polarized with 
versus without the extra d functions ((5(C1)=3.51 and 3.92, respectively). The NPA charges 
are not significantly affected ((5(C1)=2.68 in both cases), while the APT charges actually 
show a mild opposite trend (2.58 vs. 2.53). 

Summing up, the chemical significance of the extra d functions is thus that they improve 
the ability of the 3d orbital to act as an acceptor for backbonding from the oxygen lone pair 
orbitals. This situation is, in fact, somewhat reminescent of the role of the 3d orbital in the 



CaO molecule 



Ml 



As a result of all this, basis set convergence for the aug-cc-pVnZ series is quite 
unsatisfactory. Somewhat more satisfying results are obtained with the Wilson-Peterson- 



6 



Dunning aug-cc-pV(n+d)Z series, although even here, there is still an unusually large 
increment of 0.68 kcal/mol (for CI2O7: 1.29 kcal/mol) between the aug-cc-pV(5+d)Z and 
aug-cc-pV(6+d_)Z basis sets. With the aug-cc-pVnZ+2dlf series advocated in the original 
W1/W2 paper jj, convergence is rather more satisfactory: the aug-cc-pV5Z+2dlf basis set 
actually yields a larger binding energy (78.81 kcal/mol) than the aug-cc-pV(6+d)Z and aug- 
cc-pV6Z+d basis set (78.70 and 78.76 kcal/mol, respectively). Our best directly computed 
value, with the aug-cc-pV6Z+2dlf basis set, is only 0.05 kcal/mol greater. Note that the 
contribution of the high-exponent / function decays rapidly in the aug-cc-pVnZ+2dlf series: 
from 0.98 kcal/mol for n=T to 0.02 kcal/mol for n=6. Note also that even for the aug-cc- 
pV6Z+d basis set, adding a second hard d function still contributes 0.08 kcal/mol. 

Turning to Jensen's polarization consistent basis sets — which were specifically developed 
for SCF or DFT applications rather than for correlated ab initio calculations — we still see 
a 0.58 kcal/mol increment between uncontracted aug-pc3+d and aug-pc4 basis sets. (Jensen 
recommends using uncontracted basis sets for SCF energy extrapolations.) Using Jensen's 
recommended extrapolation formulaj3] from aug-pc2+2d, aug-pc3+d, and aug-pc4 data 
(see footnote c to Table 1 for details), we obtain a best estimate for TAE(SCF)=78.94 
kcal/mol: for comparison, the best directly computed values are 78.70 kcal/mol for aug-cc- 
pV(6-|-d)Z, 78.84 kcal/mol for aug-cc-pV6Z+2d, 78.86 kcal/mol for aug-cc-pV6Z+2dlf, and 
78.77 kcal/mol for uncontracted aug-pc4. 

Using the aug-cc-pV{5,6}Z+2dlf data instead with the A + B/L^ extrapolation used in 
W2 theory, we obtain 78.89 kcal/mol as our limit: taking into account Jensen's observation 
that the SCF energy for such large basis sets appears to converge much faster than (see 



also Ref. 



45j), we take the average of this extrapolation and the raw aug-cc-pV6Z+2dlf 



value, to finally obtain TAE(SCF)=78.88 kcal/mol. The Wlw extrapolation falls far short 
at 76.35 kcal/mol, while the performance of the W2w extrapolation, 78.61 kcal/mol, could 
be deemed acceptable. (With the " -|-2dlf ' basis set series of standard W2 theory, we would 
obtain 79.09 kcal/mol, which overshoots the true limit.) Using the efficient direct SCF codes 
presently available, it might well be that the preferred approach for the SCF component in 
problematic cases like this — rather than an extrapolation on purely empirical grounds — 
would be to carry out a single straight SCF calculation with the largest routinely feasible 
basis set. Even in the present extreme case, using the raw SCF/aug-cc-pV6Z+2dlf or 
SCF/aug-pc4+d numbers causes errors under 0.1 kcal/mol. 
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An SCF/aug-cc-pV6Z+2dlf calculation on CI2O7, which requires no less than 1,743 
basis functions, yields TAE(SCF)=-30.81 kcal/mol. A + B/L^ extrapolation yields -30.76 
kcal/mol, averaging as for HCIO4 -30.79 kcal/mol. (The raw SCF/aug-cc-pV(6+d)Z and 
SCF/aug-pc4 values are in error by 0.31 and 0.19 kcal/mol, respectively, suggesting aug- 
cc-pV6Z+2dlf to be the basis set of choice for "single point" SCF limit calculations.) In 
terms of convergence along basis set sequences, the same trends as for HCIO4 are seen in 
amplified form. We note in passing that the SCF reaction energy of the isodesmic reaction 
2 HCIO4 CI2O7 + H2O converges quite rapidly with the basis, as long as at least one 
hard d function is present. 



B. Valence correlation energy 

As noted previously for SO2 and SO-i.jl^. basis set convergence for the valence 
correlation component is not anomalous for molecules with severe inner polarization. As 
for any molecule with highly polar bonds, however (see, e.g., SiF4 0, BF3 Q, and the 
like), a significant difference between Wlw and W2w is to be expected. In the case of HCIO4, 
the discrepancy for the CCSD correlation energy is quite modest at 0.5 kcal/mol, and it is 
partly compensated by a discrepancy of -0.2 kcal/mol for the connected triple excitations 
contribution. 

Even a CCSD/aug-cc-pV(Q+d)Z calculation on CI2O7 was only barely feasible with the 
available hardware: CCSD(T)/aug-cc-pV(Q+d)Z and CCSD/aug-cc-pV(5-|-d)Z calculations 
are plainly impossible. Considering the good agreement between Wl and W2 for the 
valence correlation contributions, however, we can safely assume that this component is 
well reproduced for the isodesmic reaction. 

Neither molecule shows any obvious sign of severe static correlation, neither from the 
Too diagnostics^] (0.019 for IICIO4, 0.021 for CI2O7), nor from the largest coupled cluster 
amplitudes. Yet the fact that three-quarters of the atomization energy of HCIO4 results from 
correlation, and that CI2O7 is actually slightly metastable at the Hartree-Fock level, suggest 
that higher-order correlation effects could be somewhat important for these molecules. 
Unfortunately, a W3 calculation for HCIO4 — let alone CI2O7 — is absolutely impossible 
with the present state of technology. This appears to be the single greatest source of 
uncertainty in our calculations. 
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C. Inner-shell correlation energy 



Like with other second-row molecules, despite the rather large absolute inner-shell 
correlation energies, the inner-shell contribution to the molecular binding energy is rather 
small (0.90 kcal/mol for HCIO4). A CCSD(T)/MTsmall calculation on CI2O7 is not feasible. 
We were able to carry out a CCSD(T)/cc-pwCVTZ calculation, but this basis set is clearly 
woefully inadequate for HCIO4, recovering only about half the contribution. We might be 
able to rely on error compensation for the isodesmic reaction, and the inner-shell correlation 
contribution to its reaction energy is found to be -0.18 kcal/mol at the CCSD(T)/cc- 
pwCVTZ level. From the higher- level results for HCIO4 and H2O, we can then extract a 
best estimate of 1.38 kcal/mol for the inner-shell contribution to the CI2O7 binding energy. 

D. Relativistic effects 

The original Wl and W2 protocols called for mass- velocity and Darwin corrections from 
ACPF/MTsmall calculations. W3 theory, as well as the Wl and W2 implementations in 
the popular Gaussian 03 program package, use CCSD(T) calculations with the Douglas- 
Kroll approximation to obtain the scalar relativistc contribution. We have followed the 
same (more rigorous) approach here, as the scalar relativistic contributions are generally 
significant for this type of molecule (second-row atoms in high oxidation states surrounded 
by strongly electronegative elements). At the DK-CCSD(T)/MTsmall level, we find -2.69 
kcal/mol for HCIO4 and a hefty -4.80 kcal/mol for CI2O7. For HCIO4, we recalculated the 
contribution using 

AVQZ quality. We obtain fundamentally the same result, and upon extrapolation to the 
infinite-basis limit, we actually reproduce the MTsmall result to two decimal places. We 
therefore have not attempted any DK-CCSD(T)/AVQZ calculation for CI2O7: for the sake 
of completeness, the DK-CCSD(T)/AVTZ value is -4.73 kcal/mol. 

The first-order spin-orbit contribution results exclusively from the atomic fine structures 
and requires no further comment. We expect second-order spin-orbit coupling corrections 
to be on the order of 0.1 kcal/mol or less, well below the more important potential error 
sources in these calculations. 
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E. Zero-point and thermal corrections 



The scaled B3LYP/cc-pV(T+d)Z harmonic frequencies normally used for the zero-point 
energy in Wlw theory is clearly falling short a bit in these cases. From a very recent 
DFT anharmonic force field study 49] on both molecules, we have available anharmonic 
ZPVEs at the B97-l/aug-pc3+d level. Our "best estimate" value for HCIO4 was obtained 
by combining the anharmonic ZPVE with one-half the difference between the computed and 
observed fundamental frequencies. For CI2O7 the assignment of the experimental vibrational 
spectrum is too fraught with ambiguities to allow for the same approach: here we have 
assumed the ratio between the "best estimate" and directly computed ZPVEs for HCIO4 to 
be transferable to CI2O7. 

The thermal correction represents the minor complication of hindered CIO3 rotations. In 
the case of HCIO4, we find the internal rotation barrier V to be 0.647 kcal/mol at the B97- 
l/aug-pc2+2d level, which corresponds to a V/RT ratio of nearly one. At the same level, 
we obtain a reduced moment for the internal rotor of 2.9853 amu.bohr^. With a de-facto 
rotor periodicity of three, we obtain 1/Qf = 0.52 from eq.(4) of Pitzer and Gwinnj^. By 
interpolating their Table V, we obtain 1.054 e.u., or 0.53 R, for Eintrot/T, i.e., very close 
to the free rotor limit of R/2 and considerably removed from the low-frequency harmonic 
vibration limit of R. Our 'best estimate' values reflect this correction for HCIO4, and twice 
this correction for CI2O7. Clearly, it will cancel in the isodesmic reaction. 



F. Heats of formation 



Our final best estimate for the heat of formation of HCIO4 is -0.6 kcal/mol. We expect 
this value to be accurate to about 1 kcal/mol, with the largest potential source of error 
being post-CCSD(T) correlation effects. 

This value actually agrees quite well with Colussi and Grela's group additivity estimate 
of -1.5 kcal/mol. The large discrepancy between the G2 theory 52] result of Francisco jisj] . 
AffjQ =+10.8 kcal/mol and our own absolute-zero heat of formation of 2.4 kcal/mol should 
not be surprising considering the extreme basis set sensitivity pointed out here. Using 
the more up-to-date G3 theory we obtain TAEo=307.9 kcal/mol (compared to 305.4 
kcal/mol for G2), or 306.2 kcal/mol after applying the atomic spin-orbit correction (- 
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1.73 kcal/mol): G3 thus still underbinds the molecule by 7.6 kcal/mol compared to our 
best estimate (TAEo=313.8 kcal/mol). G3X theory Q, 

which was recently developed in 
an attempt to (inter alia) reduce errors for (pseudo)hypervalent compounds, recovers an 
additional 5.64 kcal/mol worth of binding energy and thus reduces the error to a respectable 
2 kcal/mol. While G3 uses an MP2/6-31G* reference geometry (which will have very 
substantial basis set incompleteness error for this type of molecules), G3X uses B3LYP/6- 
31G(2df,2p) — in particular the extra d function very significantly affects the computed 
reference geometry. In addition, G3X involves a step where a single g function is added to 
the basis set at the HF level: this does affect the energy by 2.01 kcal/mol, which appears to 
be a slight overestimate, as we find the g function contribution in the aug-pc3 basis set to 
be only 1.45 kcal/mol. 

Our best estimated heat of formation for CI2O7 (obtained from the Wlw isodesmic 
reaction energy and our best calculated data for HCIO4 and H2O) is 65.9 kcal/mol, to 
which we attach a conservative error bar of 2 kcal/mol. This is actually quite close to 
the Wagman et al. 12] estimate of Aifj 298~65.0 kcal/mol and the Riihl et al. mass 
spectrometric value of 65±4 kcal/mol. Li et al. used a variety of methods from the G2 
and G3 family, and found 76.8 kcal/mol at the G3 level, 84.9 kcal/mol at the G3(MP3) level, 
and 66.7 kcal/mol at the G2(MP2) level, which is fortuitously in much better agreement 
with our best value than the more sophisticated G3-based methods. Once again, we find 
that G3X theory puts in a much better performance: after spin-orbit correction, we obtain 
TAEo=399.47 kcal/mol, again within 2 kcal/mol of our best estimate. The G3X g function 
is found to account for 3.94 kcal/mol (once again, an overestimate compared to the 2.86 
kcal/mol from the g functions in the aug-pc3 basis set), with most of the rest once again 
being accounted for by the superior reference geometry. 

Sicre and Coboslsil obtained 93.1 kcal/mol at the B3LYP/6-311+G(3d2f) level, 78.9 
kcal/mol at the mPWlPW91 level with the same basis set, 86.2 kcal/mol at the 
G3(MP2)//B3LYP level and 79.5 kcal/mol at the G3(MP2)//B3LYP/6-311+G(3d2f) level. 
These same authors' best estimate, from isodesmic reaction schemes, is 61.5 kcal/mol, 
appreciably below our best estimate. 
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IV. CONCLUSIONS 



The heats of formation of HCIO4 and CI2O7 have been determined to chemical 
accuracy for the first time by means of Wl and W2 theory. These molecules exhibit 
particularly severe degrees of inner polarization, and as such obtaining a basis-set limit SCF 
component to the total atomization energy becomes a challenge. (Adding high-exponent 
d functions to a standard spd basis set has an effect on the order of 100 kcal/mol for 
CI2O7.) Wilson's aug-cc-pV(n-|-d)Z basis sets represent a dramatic improvement over the 
standard aug-cc-pVnZ basis sets, while the aug-cc-pVnZ-|-2dlf sequence converges still more 
rapidly. Jensen's polarization consistent basis sets still require additional high-exponent d 
functions: for smooth convergence we suggest the {aug-pcl+3d,aug-pc2+2d,aug-pc3+d,aug- 
pc4} sequence. The role of the tight d functions is shown to be an improved description 
of the CI (3d) Rydberg orbital, enhancing its ability to receive back-bonding from the 
oxygen lone pairs. In problematic cases like this (or indeed in general), a single SCF/aug- 
cc-pV6Z-|-2dlf calculation may be preferable over empirically motivated extrapolations. 
Our best estimate heats of formation are Aifj298[HC104(g)] = — 0.6±1 kcal/mol and 
Ai/j298[C'^207(g)] = 65.9±2 kcal/mol, the largest source of uncertainty being our inability 
to account for post-CCSD(T) correlation effects. While G2 and G3 theory have fairly large 
errors, G3X theory reproduces both values to within 2 kcal/mol. 
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TABLE I: Contributions to the total atomization energies and heats of formation of HCIO4 and 
CI2O7. Ah data in kcal/mol 
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(a) reaction energy of 2 HCIO4 — > CI2O7 + H2O 

(b) Best estimate 1.44 kcal/mol if isodesmic reaction assumed to be tliermoneutral w.r.t. inner shell correlation 

(c) HCIO4: From A + B/L^ and AV{5,6}Z+2dlf data: 78.89 kcal/mol. Assuming this overestimate, and taking average of 
raw value and extrapolated, best estimate is 78.88 kcal/mol. CI2O7: -30.76 and -30.79 kcal/mol, respectively, basis sets and 
Jensen's extrapolation formula E{L, ns) = Eoc + A{L + 1) cxp{~B^/nJ), with rts taken for CI: B=4. 901885, -Boo=78.94 
kcal/mol. (Nonlinear equation solved using "goalseek" feature of Excel.) 

(d) From CCSD(T)/cc-pwCVTZ calculations. With cc-pwCVTZ basis set: 0.336 kcal/mol H20, 0.581 kcal/mol CI2O7, 0.546 
HCIO4, hence => -0.175 for reaction. Hence best estimate for CI2O7 from reaction and best available data for H2O and 
HCIO4: 1.38 kcal/mol 

(e) Expt. value: -57.80 kcal/mol 

(f) From B97-l/aug-pc3+d anharmonic force fieldf^l 

(g) From anharmonic force field and ''iCoxpt.) — i^i(calc.)=0.20 kcal/mol 

(h) From anharmonic force field I49I : ZPVE=19.30 kcal/mol. Expt. spectral assignment problematic 4^: assuming similar 
relative error as for HCIO4 yields scaling factor of 1.0115, hence best estim. ZPVE=19.52 kcal/mol 

(i) extrapolated from -2.64 and -2.67 kcal/mol, respectively, with relativistic correlation consistent basis sets (see Ref.Q) 
(j) from Gaussian calc. with hindered rotor stuff: reduced moment=3.0142 amu.bohr^. In these units, eq. (4) of Pitzer & 



Gwinn.5l| becomes Qf = 2.815\/0.004657redT/n, with rotor periodicity n de facto =3. Hence 1/Qy=0.52; with V/RT fsi 1, 
from Table V of Pitzer & Gwinn we obtain by interpolation -Eint.rot. = 0.53-R. Thus the RRHO H298 — ^^0 has to be reduced 
by 0.47 RT=0.28 kcal/mol at 298 K. 

(k) from Wl reaction energy of (a) and best calculated data for HCIO4 and H2O 
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TABLE II: Basis set convergence of the SCF contribution to the total atomization energies 
(kcal/mol) of HCIO4 and CI2O7 
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FIG. 1: B3LYP/cc-pV(T+d)Z (in square brackets) and B97-l/aug-pc3+d geometries 
and Cl207(A, degrees) 



49| for HCIO4 
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